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PURPOSE 


The  purpose  of  this  research  is  to  furnish  knowledge  that  can  be 
used  in  the  development  of  objective  methods  of  analyzing  and  fore¬ 
casting  tropical  weather  and  in  understanding  the  dynamics  of  tropical 
weather  phenomena.  It  is  planned  that  the  objective  techniques  will 
be  designed  for  electronic  computation  in  order  to  gain  speed  and 
accuracy,  and  to  reduce  personnel  requirements  in  an  operational  situation. 
The  investigations  are  divided  into  the  following  tasks: 

(1)  Analyze  the  surface  weather,  the  three-dimensional 
structure,  the  changes  with  time,  and  the  movement, 
of  selected  cases  of  representative  meteorological 
phenomena,  by  using  conventional  methods.  If 
practicable,  Tiros  data  will  be  included  in  the 
analysis . 

(2)  Investigate  and,  insofar  as  possible,  develop 
objective  analysis  techniques  applicable  to  tropical 
phenomena  and  compare  the  objective  analyses  with 
the  analyses  of  Task  (1)  above. 

(3)  Utilizing  the  analyses  of  Tasks  (1)  and  (2)  above, 
carry  out  dynamical  studies  of  such  topics  as 

the  forces  predominant  in  various  phenomena,  the 
conservation  of  fields  of  vorticity,  divergence 
and  deformation,  and  the  importance  of  orography 
and  low-level  energy  inputs. 
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ABSTRACT 


During  the  second  quarter,  IBM  card  data  for  the  Caribbean  for  the 
period  5-8  May  1959  were  processed  in  the  form  of  layer  averages  of  wind 
components,  height,  temperature,  and  relative  humidity.  The  averaged 
quantities  were  mapped  by  the  computer,  and  subjective  analyses  were 
partially  completed.  A  computer  program  for  obtaining  a  streamf unction 
was  completed.  The  streamf unction  is  made  to  fit  observed  winds,  while 
observed  heights  are  used  as  an  initial  guess.  No  boundary  conditions 
are  required.  A  preliminary  analysis  indicates  that  the  results  are 
quite  satisfactory.  Analyses  of  computed  values  of  vorticity,  divergence, 
and  vertical  motion  were  also  begun,  and  average  magnitudes  of  these 
quantities  were  determined. 
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I  PUBLICATIONS,  LECTURES,  REPORTS,  AND  C0NPIRKNCE8 


Publications,  Lectures,  and  Reports 

On  14  September  1962,  a  paper  entitled  "Some  Applications  of 
Computers  in  Tropical  Analysis"  was  presented  by  R.  M.  Bndlich  at  the 
Fourth  Conference  on  Applied  Meteorology  at  Fort  Monroe,  Virginia. 

This  paper  summarized  the  approach  which  is  being  pursued  in  the  contractual 
work  and  followed  the  material  presented  in  the  First  Quarterly  Progress 
Report . 

Conferences 

The  Chief  Signal  Officer,  Major  General  E.  Cook,  visited  the 
Institute  on  5  September  1962.  During  his  visit,  Mr.  Endlich  and  Dr. 

Ligda  outlined  the  contractual  work  and  discussed  tropical  meteorology 
with  him. 

On  12  September  1962,  a  discussion  and  review  of  the  contractual 
work  took  place  at  Fort  Monroe  between  Mr.  Endlich  of  SRI,  and  Mr. 

Raymond  Bellucci  and  Mr.  Joseph  Valsh  of  the  Atmospheric  Physics  Branch, 
USAERDL. 
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II  FACTUAL  DATA 


A.  Task  1 

The  IBM  card  data  for  the  Caribbean  for  the  period  5-8  July  1959 
were  processed  in  the  form  of  layer  averages  of  height,  temperature, 
relative  humidity,  and  wind  components.  Analyses  of  the  various  quantities 
(drawing  of  isolines)  are  partially  completed.  Examples  of  layer- 
averaged  winds  and  heights  are  shown  in  Figs.  1  and  2,  and  will  be  dis¬ 
cussed  in  more  detail  in  Sec.  II-B.  A  point  of  considerable  interest 
is  that  the  layer-averaged  humidity  fields  (see  Fig.  3)  consist  of 
fairly  large,  regular  patterns.  It  is  believed  that  this  fact  will  be 
important  in  analysis  of  clouds  and  precipitation  and  possibly  also  in 
estimating  vertical  motions  from  wind  and  humidity  data.  Since  the 
analyses  for  the  entire  series  are  not  complete,  further  discussion 
will  be  deferred  to  the  next  quarterly  report. 

B.  Task  2 

In  regard  to  the  development  of  objective  analysis  techniques,  a 
number  of  basic  computer  programs  for  the  Caribbean  have  been  completed. 

The  programs  are  indicated  in  Fig.  4.  It  is  believed  that  they  provide 
a  means  of  rapidly  obtaining  values  of  most  quantities  of  meteorological 
interest.  Output  can  be  in  tabular  form  or  printed  at  station  locations 
(i.e.,  mapped). 

In  the  First  Quarterly  Status  Report,  the  need  for  an  objective, 
numerical  streamfunction  to  replace  conventional  subjective  stream¬ 
lines  was  discussed.  During  the  present  quarter,  a  machine  program  for 
computing  a  numerical  streamfunction  was  completed.  The  results  to 
date  are  very  satisfactory.  Since  we  anticipate  that  the  technique 
will  be  described  fully  in  a  future  scientific  report,  only  a  summary 
will  be  included  at  the  present  time.  The  horizontal  wind  vector  V  is 
represented  by  a  streamfunction  $  as  follows: 

W  =  (g/f)  <*  x  V$>  (1) 
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where  $  has  the  units  of  meters.^  It  is  assumed  that  divergent  wind 
components  can  be  suppressed  by  an  areal  averaging  process.  The  equations 
relating  wind  components  and  the  streamfunction  are 

u  =  -  (g/f)  (d$/dy),  v  =  (g/f)  (d$/dx)  .  (2) 

Since  wind  observations  are  the  primary  data  used  in  tropical  analysis, 
observed  values  of  u  and  v  are  used  to  determine  values  of  d$/dy  and 
d\|r/dx,  and  thence  values  of  Consider  a  pair  of  stations  separated 
by  distance  increments  Bx  and  By.  The  streamfunction  difference  be¬ 
tween  them  is 

6$  =  (df/dx)  Bx  +  (d$/dy)  By  =  (f/g)  (vBx  -  uBy)  (3) 


where  the  appropriate  values  of  u,  v,  and  f  are  the  average  values  along 

a  line  joining  the  stations.  These  averages  may  be  taken  as  (u  +  u. )/2, 

o  i 

(v  +  v. )/2  and  (f  +  f.)/2,  where  the  subscript  o  denotes  values  at 
o  i  o  i 

the  station  of  interest  and  i  denotes  values  at  a  neighboring  station. 
Then  Eq.  (3)  becomes 


I  =1 

yo 


(fo  +  V 

2g 


<vo  +  vA)  (xQ  -  x±)  <uo  +  u±>  (yQ  -  yj 


(4) 


An  equation  of  this  form  may  be  obtained  between  station  o  and  each  of 
its  n  neighbors — i.e.,  i  =  1,  2,  . ..,  n. 


A  reasonable  estimate  of  ^  might  be  taken  as  the  average  of  the 
individual  estimates--i.e., 

n  n 


K  =  +  +  V  [<Vo  +  V  Uo  "  Xi)  ' 

(uo  +  ui}  (yo  "  yi)] 


i=l 


i=l 


(5) 


Since  f  approaches  zero  as  the  equator  is  approached,  it  is  desirable 
to  avoid  the  use  of  f  in  equatorial  regions.  This  can  be  done  simply 
by  setting  W  1  Ik  x  vt  (where  y  has  the  units  m^sec”1')  or  by  placing  a 
fictitious  lower  limit  on  f. 
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However,  experience  in  objective  analysis  in  mid- latitudes  has  shown 

that  observations  should  be  weighted  inversely  with  their  distance  from 

the  point  of  interest. 2 Therefore,  one  estimate  ($  ),  is  obtained 

o  l 

by  considering  all  stations  lying  within  a  radius  of  5  degrees  of 

latitude  and  a  second  estimate  ($  )„  is  obtained  by  considering  stations 

o  < 

lying  between  5  and  9  degrees  of  latitude  away.  The  two  estimates  are 
then  combined  as 


*o  =  <l)l  +  ^2  <i>2  (6) 

where  |i  and  |i  are  weighting  factors  presently  set  at  0.8  and  0.2.  The 

entire  field  of  observing  stations  is  scanned  in  this  manner,  producing 
* 

adjusted  values  of  ijr  which  are  immediately  substituted  for  the  previous 

ones  according  to  the  Liebmann  iteration  process.  The  initial  values 

of  |  at  each  station  are  taken  equal  to  the  reported  heights  Z  .  Several 
1  1  * 
scans  of  the  field  are  performed  until  nearly  convergent  values  of  \|r  are 

obtained.  The  average  of  the  computed  values  of  streamf unction  remains 

approximately  equal  to  the  average  of  the  heights  used  as  the  initial 

guess.  It  was  assumed  at  the  outset  that  the  divergent  wind  components 

* 

would  be  suppressed.  Since  the  computation  of  \jr  at  each  station  depends 
upon  several  winds  out  to  a  radius  of  9°  of  latitude,  the  divergent 
features  are  automatically  suppressed  by  the  computational  procedure. 

The  direct  streamfunction  computation  described  above  appears  to 
have  several  advantages  over  other  possible  formulations.  The  general 
approach  is  simple  in  principle.  The  use  of  boundary  conditions  is 
avoided.  It  is  not  necessary  to  use  a  rectangular  grid.  The  computed 
values  of  f  can  be  compared  directly  with  observed  heights.  Wind  estimates 
made  from  satellite  data  could  be  utilized  in  the  system.  At  present, 
total  machine  time  for  computing  $  in  each  of  six  layers  and  printing 
out  the  computed  values  in  mapped  form  is  about  10  minutes  using  the 
Burroughs  220  computer. 


^References  are  listed  at  the  end  of  the  report. 
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A  qualitative  assessment  of  the  computed  streamfunction  can  be 

made  by  inspection  of  Figs.  1  and  2,  which  show  layered-averaged  winds, 

heights,  and  streamfunction.  Layers  II  and  V  are  representative  of  the 

lower  and  upper  troposphere,  respectively.  Note  the  agreement  of  winds 

and  streamfunction  in  contrast  to  the  lack  of  reasonable  relationships 

between  winds  and  height  contours.  (Twelve  hours  after  the  charts  shown, 

the  correspondence  of  heights  and  winds  was  even  less  satisfactory.)  A 

preliminary  quantitative  evaluation  can  be  made  from  the  information 

contained  in  Table  I,  The  first  quantity  is  the  root-mean-square  (RMS) 

value  of  the  observational  error  in  height  (Z) .  The  second  quantity  is 

* 

the  RMS  value  of  the  difference  between  Z  and  the  streamfunction  i|r.  The 
other  four  quantities  are  RMS  values  of  12-hour  and  24-hour  changes  of 
height  and  streamfunction.  It  is  interesting  to  note  that  the  RMS  12- 
hour  height  changes  exceed  the  24-hour  changes,  contrary  to  expectations. 

Table  I 

COMPARISON  OF  HEIGHT  (Z)  AND  STREAMFUNCTION  (1) 

DURING  THE  PERIOD  00  GOT,  5  MAY  1959  TO  00  GOT,  6  MAY  1959 
IN  THE  CARIBBEAN  (30  STATIONS) 


Layer  II  RMS  Error  of  Z  »  5  m 

RMS  (Z  -  |)  =  9  m 

RMS  12-hour  height  change  =  13  m 

RMS  24-hour  height  change  =  9m 

RMS  12-hour  \|f  change  =  5  m 

RMS  24-hour  f  change  =  9m 

Layer  V  RMS  Error  of  Z  ~  20  m 

RMS  (Z  -  |)  =  27  m 

RMS  12-hour  height  change  s  31  m 

RMS  24-hour  height  change  =  26  m 

RMS  12 -hour  $  change  :  15  m 

* 

RMS  24-hour  \jr  change  =  17  m 
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The  reason  for  this  result  has  not  yet  been  investigated.  However,  the 
main  point  of  interest  is  a  comparison  of  streamfunction  changes  and 
height  changes.  The  smaller  values  of  streamfunction  changes  indicate 
that  streamfunction  is  a  considerably  more  persistent  quantity  than 
height.  In  fact,  in  Layer  V  the  RMS  streamfunction  change  is  slightly 
smaller  than  the  RMS  error  of  height.  These  results  are  very  encouraging. 

As  indicated  in  Fig.  4,  values  of  vorticity,  divergence,  deformation, 
and  vertical  motion  have  been  computed  from  wind  data  by  use  of  the  tri¬ 
angle  technique.3  Examples  of  the  computed  values  in  Layers  II  and  V 
are  shown  in  Figs.  5  and  6  for  purposes  of  illustration.  In  Layer  II, 
the  relative  vorticity  (solid  lines)  is  negative  in  the  anticyclone  (Fig. 
1)  and  positive  in  two  centers  related  to  the  trough  in  the  southeastern 
portion  of  the  Caribbean.  The  divergence  field  contains  four  centers 
of  size  sufficient  to  insure  that  they  are  real.  The  vertical  motion  at 
the  top  of  Layer  II  is  shown  as  part  of  Fig.  3.  The  vertical  motion  is 
predominantly  negative  (downward)  except  in  the  region  of  the  trough, 
and  appears  to  have  a  reasonable  correspondence  with  the  humidity  field. 

In  Layer  V,  the  vorticity  field  (Fig.  6)  consists  of  four  main 
centers  which  appear  to  characterize  important  features  of  the  flow 
shown  in  Fig.  2.  The  northwesterly  jet  lying  west  of  Florida  has  cyclonic 
relative  vorticity  on  the  low-pressure  side  and  anticyclonic  vorticity 
on  the  other  side.  The  streamline  trough  is  a  region  of  cyclonic  vor¬ 
ticity,  and  the  south  side  of  the  jet  lying  east  of  the  trough  is  a 
region  of  anticyclonic  vorticity.  The  divergence  field  in  Layer  V  con¬ 
tains  a  major  center  of  convergence  over  the  Gulf  of  Mexico  in  advance 
of  an  isotach  maximum.  The  other  centers  are  apparently  related  to  the 
trough  and  jet  stream. 

At  this  time  we  have  not  attempted  to  evaluate  the  computed  values 
of  divergence,  vorticity,  etc.,  nor  to  determine  their  proper  use  in 
analysis  and  forecasting.  These  subjects  require  a  considerable  amount 
of  further  study. 
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FIG.  5  RELATIVE  VORTICITY  AND  DIVERGENCE  IN  LAYER  II 


c. 


Task  3 


The  average  magnitudes  of  several  terms  important  in  dynamical 
equations  have  been  computed  from  data  for  5-6  May  1959  and  28  July 
1959.  Discussions  at  the  Conference  on  Tropical  Meteorology4  revealed 
differing  opinions  regarding  these  magnitudes.  Table  II  lists  the  re¬ 
sults  of  our  computations  for  the  Caribbean. 

Table  II 

AVERAGE  MAGNITUDES  OF  DYNAMICAL  TERMS 


Divergence  0 . 9 
Vorticity  1.1 
Resultant  Deformation  1.4 


Layer  II 

Layer 

-5  -1 

x  10  Bee 

1.5 

x  10-5 

x  10"5 

3.0 

x  10-5 

x  10~5 

3.2 

x  10-5 

It  is  rather  surprising  that  these  values  are  nearly  as  large  as  typical 
values  computed  for  mid- latitudes . 3  This  may  be  due  in  part  to  the 
presence  of  unusually  active  circulations  in  the  Caribbean  during  the 
periods  under  investigation. 

Consideration  has  been  given  to  several  other  topics,  including 

the  determination  of  synoptic-scale  vertical  motions  from  relative 

humidity  changes  experienced  by  parcels,  objective  techniques  of  kine- 

matical  forecasting,  methods  of  applying  the  vorticity  equation  in  the 

tropics,  and  use  of  the  conventional  balance  equation  to  determine  Z 
* 

from  ijr.  These  topics  will  be  discussed  in  later  reports. 
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Ill  PROGRAM  FOR  THE  NEXT  QUARTER 


The  analyses  for  the  period  5-8  May  1959  will  be  completed  and  a 
major  effort  will  be  expended  in  evaluation  of  the  objectively  computed 
streamf unction,  divergence,  vertical  motions,  etc.  The  effects  of 
scale  (areas  considered)  upon  computed  values  will  be  investigated. 
Further  studies  related  to  dynamical  quantities  will  be  carried  out. 

It  is  also  expected  that  a  scientific  report  describing  the  stream- 
function  computation  will  be  begun.  At  this  time,  no  definite  plans 
for  travel  during  the  quarter  have  been  made. 
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Approximately  55  hours 
Approximately  20  hours 
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